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Abstract 

We consider the classical dynamics of a particle in a one-dimensional space- 
periodic potential U{X) = U{X + 27r) under the influence of a time-periodic 
space- homogeneous external field E{t) = E(t + T). If E{t) is neither sym- 
metric function oft nor antisymmetric under time shifts E{tzizT/2) ^ —E{t), 
an ensemble of trajectories with zero current at t = yields a nonzero finite 
current as t — > oo. We explain this effect using symmetry considerations and 
perturbation theory. Finally we add dissipation (friction) and demonstrate 
that the resulting set of attractors keeps the broken symmetry property in 

the basins of attraction and leads to directed currents as well. 
05.45.-a, 05.60.Cd, 05.45.Ac 
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I. INTRODUCTION 



Transport phenomena are at the heart of many problems in physics. Nonhnear effects 
(as well as their quantized counterparts) may lead to many novel results in this area even 
for seemingly simple models (see e.g. [|l|). Well-known applications include the dynamics of 
Josephson junctions and electronic transport through superlattices 0, to name a few. 
In the bulk of theoretical work on transport phenomena nonzero dc currents are obtained 
by applying time-dependent fields with nonzero mean. It is normally expected that the 
opposite case may not lead to a nonzero dc current. However it has been also known for 
a long time that nonlinear dynamical systems may allow for generation of ac fields from 
external dc fields (Josephson effect) and even vice versa 0. Of course what matters is a 
proper average over initial conditions, so that one has to ask whether there exist (or do not 
exist) sets of solutions which cancel their contribution to the total current. This question 
calls for an analysis of the symmetry properties of the system under consideration. 

Let us make things more precise by considering a paradigmatic equation of the following 
type: 

X + -fX + f{X) + E{t) = . (1) 

Functions / and E are bounded and periodic with period 27r and T = 2tt/uj respectively 
and have zero mean, and max(|/(X)|) ~ 1. This equation describes e.g. a particle moving 
in a periodic potential U{X) with f{X) = U'{X) in one space dimension under the influ- 
ence of a periodic external field with friction 0. It also may describe the current- voltage 
properties of a small Josephson junction under the action of a time-periodic external current 
(here X becomes the phase difference of the complex order parameter across the junction). 
This equation has been considered by numerous authors, however typically with harmonic 
functions / and E. We will show below that this choice induces symmetries which lead to 
zero total dc current. The purpose of this letter is to demonstrate that a proper lowering of 
the symmetries of even E{t) alone (still keeping its above defined properties) will lead to a 
nonzero dc current. 
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II. DISSIPATIONLESS CASE 7 = 



We first consider the case of zero friction 7 = in (^. In the hmit of large velocities 
|X| > 1, /(X) can be neglected and the solution X{t) = Xq + Pot + Jodt' J^' E{t")dt" has 
a bounded first derivative. Thus the time average over the velocity on a given trajectory is 
a well defined nondiverging quantity. 

To characterize the relevant symmetries of (0) we have to consider transformations in 
X, t which lead to a change of sign in P. These are i) a reflection in X — —X and a 
shift in t or ii) a shift in X and a reflection in t ^ —t. We need first to characterize the 
relevant symmetries of /(X) and E{t). For that we expand / and E into Fourier series: 
f{X) = Y.kfk^^^'^ ! E{t) = J2kEke"^''^- Zero mean implies /o = -Eq = 0, and reality yields 
fk = = E*_^, {A* means complex conjugation). If /(X) = U\X) is antisymmetric 

after some appropriate argument shift f{X + X) = —f{—X + X) we call /(X) possessing 
symmetry. If E{t) is symmetric after some appropriate argument shift E{t + T) = E{—t + T) 
we call E{t) possessing Eg symmetry. If E{t) changes sign after a fixed argument shift (which 
trivially can be only equal to any odd multiple of T/2) E(t) = —E{t + T/2), resulting in 
E2k = 0, we call E(t) possessing Egh symmetry. 

Now we can define the two relevant symmetry cases of (|1]) called Sa and Sb below. If 
functions /(X) and E{t) possess fa and Egh symmetries respectively, then (|I|) is invariant 
under symmetry Sa- X {—X + 2X), t — > t + T/2. If function E(t) possesses Eg symmetry, 
(H) is invariant under symmetry Sb'. t — > (— t + 2r). 

Given a trajectory X(t;Xo,Po)) Pit'i^o^Po) with X(to;Xo,Po) = -'^o and 
P{to; Xq, Pq) = Pq the presence of any of the two symmetries Sa, Sb allows to gener- 
ate new trajectories given by 

Sa-. -X{t + T/2;Xo,Po) + 2X , -P (t + T/2; Xq, Pq) , (2) 
Sb-. X(-t + 2r;Xo,Po) , -P(-t + 2r; Xq, Pq) . (3) 

Note that these transformations change the sign of the velocity P. Consequently the time 
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average of P on the original trajectory will be opposite to the time averages of P on the 
generated new trajectories. There can be more symmetry operations generating other tra- 
jectories, but those will not change the sign of P and are thus not of interest here. 

The dynamical evolution of (|1|) allows both for quasiperiodic solutions (cyclic in X for 
large Pq and periodic in X for small Pq) and chaotic trajectories embedded in a stochastic 
layer [|l|]. Assuming that ergodicity holds in the stochastic layer we conclude that the average 
velocity will be one and the same for all trajectories of the layer. Since Sa and Sb when 
applied to a trajectory inside the layer generate again trajectories inside the layer, the 
presence of any of these symmetries implies that the time-averaged velocity of any trajectory 
in the layer will be zero. Note that we cannot obtain such a conclusion if both symmetries 
are absent! Indeed in Fig. 1 we show the long-time run X{t) for a trajectory in the layer for 
several cases with and without symmetries Sa, Sh- While with Sa, Sh we find zero average 
velocities, we observe that the loss of Sa, Sb leads to a nonzero average velocity which is 
independent on the initial conditions but whose sign depends on the way the symmetry 
is broken. The dynamics is characterized by anomalous transport, i.e. by Levy flights of 
different length interrupted by direction-changing perturbations. Nonzero current appears 
due to a desymmetrization between Levy flights to the left and right, respectively. Especially 
trajectory 2 in Fig. 1 yields a nonzero velocity for a spatially symmetric U{X). 

To answer the question of how to invert the direction of a nonzero current in the stochastic 
layer, we note that considering the equation X + f{X) +E{—t) = we arrive back at (|l]) by 
substitution t' = —t. So the current can be inverted by applying E{—t) instead of E{t) in 
(0). A second way is to consider equation X — f{—X) — E{t) = which after substitution 
X' = —X again is mapped onto (|l|). Thus another way of inverting the current is to apply 
—f{—X) instead of f{X) and —E{t) instead of E{t) in (|I|). There is no simple way to invert 
the current by just inverting space i.e. by considering f{—X). 

To get a grasp of this result we consider the quasiperiodic cyclic regime for U{X) = 
— cosX and E{t) = Ei cosujt + E2 cos{2ujt + a). Note that Sa symmetry is present if £'2 = 
OT El = and Sb symmetry is present if a = 0, vr or i?i = or i?2 = 0. Each individual 
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trajectory for sufficiently large Pq gives a nonzero average velocity. The question is whether 
we obtain a nonzero velocity after averaging over initial conditions with some distribution 
function p(Xo,Po5^o) reflecting equilibrium properties, at least of course p(Xo, Pq, ^o) = 
p{Xq, —Pojto). Here to is the time when the trajectories with initial conditions Xq,Pq are 
started. In the simplest case we might assume that p is independent of to- Consider the 
case Pq ^ I and u ^ Pq. In that case we can separate the solution X{t) into a slow part 
Xs{t) and a small fast part ^(t). Expanding to linear order in the fast variable yields 

Xs + ^- sin X, - cos(X,)e + E{t) = . (4) 

Collecting the fast variables we find ^ — cos{Xs)^ + E(t) = . This equation has to be solved 
by assuming that Xg is constant and skipping the slow homogeneous solution part. We find 
^ = Ai cos ujt + A2 cos{2ujt + a) with Ai = —Ei/[uj'^ — cosXs] and A2 = —E2/[4:Uj'^ — cosXs\. 
Final averaging over the fast variables in gives Xg — sinX^ = 0. The crucial point is to 
observe that the initial condition is now Xq = Xs(to)+^(to) , Pq = -'^s(to)+^(^o)- Since ^{t) 
is a completely defined function, defining the initial conditions for X, P we obtain initial 
conditions for the slow variables. The symmetry breaking will be hidden there. Indeed, 
averaging over time we find < P{t) >=< Xgit) >. Assuming e.g. large values of Pq the 



time average velocity of the slow variable will be simply < Xs{t) >= sgn{Po)^/2H^[l — 
l/iAHl) + 0(Po~^)] with 2H, = P^ - 2cosX,. Expanding < X,(t) > in powers of I/Pq 
we will encounter terms PQ"^^^(to) cos^[Xo — ^(to)]. Averaging over Xq and to we obtain in 
leading order for the average velocity 

which remains nonzero and will contribute to an average nonzero current after further av- 
eraging over Pq. Note that the directed current disappears if _Ei = or P2 = or a = , vr 
when the mentioned symmetries are restored. The current direction is defined in this per- 
turbation limit by the sign of the product E2 sin a. Finally in the limit Po ^ 00 the current 
amplitude tends to zero, although the symmetries are not restored. The reason is that in 



this limit we recover the problem of a free particle moving under the influence of an external 
field E{t) which can be easily solved |0. Averaging over to in this case yields zero total 
current. It follows that nonzero total currents occur if symmetries Sa and Sb are violated 
and if we provide a mechanism of mixing of different harmonics as it happens in nonlinear 
equations of motion (see also [Q). 

We checked the above statements of the perturbation theory for the quasiperiodic regime 
by computing numerically the average velocity < > for two initial conditions with op- 
posite initial velocities iPo, taking their half sum, and finally averaging over all possible 
initial positions Xq and over the initial time to. We observe a nonzero current except for the 
symmetric values of a. Finally we did the same direct computation in the initial equation 
(HI). The results are similar. 

In order to keep the dc current nonzero the value of a should be kept fixed with time, or 
at least to be allowed to fluctuate only with small amplitude. Additional averaging over a 
will lead to a disappearance of the dc current. To our understanding this should not pose a 
technical difficulty, since one can take a monochromatic field source, and then experimentally 
generate a second harmonic from it such that the phase a is fixed. 

III. THE CASE WITH DISSIPATION 

Consider now a small but nonzero value of 7 in ([^)(see 0). Generically the phase space 
of the system will separate into basins of attraction of low-dimensional attractors. There 
exist strong hints that when being close to the Hamiltonian case, these attractors will be 
periodic orbits or limit cycles (cyclic in X) 0. The stochastic layer is transformed into 
a complex transient part in phase space, where the basins of attraction of different limit 
cycles are entangled in a complicated way. For stronger deviations from the conservative 
limit the periodic attractors undergo (period doubling) bifurcations, and finally possibly 
chaotic attractors are generated, which are however not directly related to the stochastic 
layer of the conservative limit (see also |I|). 
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Of the two symmetries Sa, Sb in the conservative case only Sa may survive for nonzero 
dissipation. Consider such a case when (H) holds. Suppose we find a limit cycle which is 
characterized by X{t + T) = X{t) + 27rm and P{t + T) = P{t), m & Z. Due to the external 
time-periodic field E{t) we have T = n2'K/uj, n E Z. The average velocity < P >= 7^ Xdt 
on such a cycle will be given by < P >= uom/n. Due to the required symmetry there will 
be also a limit cycle with < P >= —uom/n. Moreover the symmetry presence also implies 
that the basins of attraction of the two symmetry related limit cycles are equivalent. 

Assume now that we violate Sa- The two cycles previously related by symmetry to each 
other will generically continue to exist, but there is no obvious symmetry which relates them 
to each other. However after computing the average velocities, we will still find that they 
equal each other up to a sign! The symmetry breaking is in fact hidden in a desymmetrization 
of the two basins of attraction. It is this asymmetry which after averaging over initial 
condition distributions (symmetric in P) will lead to a different number of particles attracted 
to both cycles and thus to a nonzero current. To observe the desymmetrization of the basins 
locally we may tune some parameter of the equation to such a value that one of the cycles 
becomes unstable. In that case its basin of attraction shrinks to zero and disappears. If the 
other (previously symmetry related) cycle will be still stable, i.e. if its basin of attraction still 
exists, the asymmetry in the basins becomes obvious - one of them completely disappeared, 
the other one still exists. We tested these predictions and found complete agreement. We 
used 

/(X) = sinX + i;2sin(2X + 0.4), (6) 
E{t) = El sin ut + E2 sm{2ujt + 0.7) (7) 

with 7 = 0.005 and u = 1.1. The two symmetry related limit cycles {n = 1 and m = ±1) 
have been computed with a Newton method (see e.g. |P) for V2 = E2 = 0, Ei = —2.0. Then 
the parameters were changed to V2 = 0.02, Ei = —2.017 and E2 = —0.06051 and the two 
limit cycles were traced again with a Newton method. Finally the eigenvalue problem (3x3 
matrix) of the linearized phase space flow around each of the cycles has been evaluated in 
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order to check the stabihty (see for details). For the given parameter values the m = — 1 
cycle is stable (all Floquet eigenvalues inside the unit circle) while the m = 1 cycle is unstable 
(one Floquet eigenvalue is outside the unit circle). 

To observe the effect of asymmetry of basins of attraction globally, we computed the 
ensemble averaged velocity for a distribution of initial conditions in the phase space of (|ip 
with forces (^-(0). The distribution was uniform in X and to (40 points on the interval from 
to 27T for each of them) and 2 x 20 points symmetrically chosen on the P-axis according 
to a Maxwell distribution with inverse dimensionless temperature f3 = 0.01. In total 64000 
trajectories have been computed. The velocity per trajectory averaged over the whole set of 
trajectories is shown in Fig. 2 as a function of time for the case with Sa symmetry (curve 1) 
and the one without Sa symmetry (curve 2). While the first case gives zero current density 
as if: — s> oo, the second case yields nonzero negative current density in this limit. 

In order to invert the direction of a nonzero total current we have to apply — /(— X) 
instead of f{X) and —E{t) instead of E{t) in (|I|). In contrast to the dissipationless case 
we cannot just invert time in E{t) but have to perform a combined transformation both 
in space and time. Taking just f{—X) or E{—t) may or may not lead to a change of the 
current direction. Recall that directed currents can be generated by keeping U (X) = U{—X) 
and lowering the symmetry in E(t) only. In that case the current direction is inverted by 
applying —E(t). 



IV. DISCUSSION 

There exist a lot of publications on the properties of (|l|) with 7 = (and similar equations 
reduced to discrete maps), however we did not find studies of such a system when both 
symmetries Sa and Sb are broken. Evidently, when taking / and E with only one harmonic, 
no symmetry broken transport is possible. The closest study in this respect we found in 
p!0| , where however, as explicitely stated, the symmetry was kept, leading to zero current 



when averaging over all possible trajectories. The overdamped case was studied in Wl . 
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Finally we want to discuss the relation of our results to the well-known case of directed 
currents for particles moving in so-called ratchet potentials under the influence of friction and 
a stochastic force (see [jl2| and references therein). These potentials lack inversion symmetry 
in space and thus lack fa symmetry (see above). However the noise process characterizing 
the stochastic force has to be non- white (see |TB[ for details) . It was then found that proper 
correlations in the noise allow for directed currents even in the presence of fa symmetry, 



i.e. for "non-ratchet" potentials. In [|14| these equations have been modified by adding 
time-periodic fields. Note that our model allows for an easy treatment of the symmetry 
analysis, since the symmetry breaking is not hidden in higher order moments of distribution 
functions. 

If we consider corresponding quantum systems, the symmetry breaking will be reflected in 
the properties of the eigenstates, and nonzero currents can be expected as well. The addition 
of e.g. particle-particle interaction or noise can only affect the amplitude of the current, 
since the broken symmetries cannot be restored by additional interactions. Applications 
of similar ideas to coherent photocurrents in semiconductors have been reported in |T3|,|TB[. 
Further applications may include driven Josephson junctions or superlattices, electrons in 
time-dependent magnetic fields to name a few. Note that it should be much easier to realize 
experimentally our proposed symmetry breaking rather than to prepare correlated noise as 
proposed for ratchet transport. 
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FIGURE CAPTIONS 



Fig.l 

Dependence X(t) versus t for different realizations of (j^) and 7 = with f{X) = cosX + 
t;2Cos(2X + 0.4), E{t) = Ei sin(tut) + Eg sin(2tut + 0.7) and u = 2.4. (1): V2 = , Ei = 
-2.4 ,^2 = 0; (2): V2 = , E^ = -2.4 , E2 = -1.38; (3): = 0.6 , E^ = -2.4 , E^ = 
— 1.38; (4): same as (2) but with f{—X) instead of f{X). Note that in this case the direction 
of the current is not inverted as explained in the text. 

Fig.2 

The averaged velocity (see text) as a function of time for (|1]) with 7 = 0.1, u; = 2.4 and 
El = —5.23. (1): symmetric case, f 2 = £'2 = 0; (2): asymmetric case, V2 = 0.6, E2 = —5.23. 
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